Unified treatment of fluorescence and Raman scattering processes near metal surfaces 
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We present a general model study of surface-enhanced resonant Raman scattering and fluorescence 
focusing on the interplay between electromagnetic effects and the molecular dynamics. Our model 
molecule is placed close to two Ag nanoparticles, and has two electronic levels. A Franck-Condon 
mechanism provides electron-vibration coupling. Using realistic parameter values for the molecule 
we find that an electromagnetic enhancement by 10 orders of magnitude can yield Raman cross- 
sections an of the order 10 -14 cm 2 . We also discuss the dependence of an on incident laser intensity. 



PACS numbers: 33.20.Fb, 33.50.-j, 42.50.-p 

Discovered nearly three decades ago, surface-enhanced 
Raman scattering (SERS) has developed into an ex- 
tremely sensitive spectroscopic technique with, in some 
cases, single molecule sensitivity 0, 0, IE 0- It is wen 
known that SERS, as well as a range of related surface- 
enhanced optical processes, mainly results from electro- 
mag netic (EM) effects HE|. 

Although the Raman scat- 
tering cross-section <jr for a molecule in free space is 
very small (of the order of 10~ 30 cm 2 for nonresonant, 
10~ 24 cm 2 for resonant scattering), the same molecule 
placed between two metal particles may well have an 
effective ctr that is 10-12 orders of magnitude larger. 
The reason is that EM fields are strongly modified near 
and, in particular, between, metallic particles, so the lo- 
cal excitation field induced by an incident wave is much 
stronger there (by a factor that we denote M) than in 
free space. Likewise, by virtue of electromagnetic reci- 
procity, the amplitude of the radiation sent out from a 
source near the particles is equally enhanced compared 
with a source in free space. Consequently, quantities such 
as the absorption cross-section for a molecule increase by 
a factor |M| 2 , whereas Raman scattering, which involves 
both an absorption and an emission event, increases by 
a factor |A/| 4 . 

A large number of theory papers on surface- 
enhancement phenomena have focused on the electro- 
magnetic aspects, but only a few have considered the 
molecular dynamics in more detail, see, for example, 
Ref. In this work we present a general model that 
treats the electromagnetic and molecular aspects on an 
equal footing. It includes photon-molecule coupling, cou- 
pling between electronic and vibrational degrees of free- 
dom on the molecule, and radiative and non-radiative 
damping mechanisms 8] , and is analyzed by means of 
a density-matrix calculation. The model lets us study 
not only how the molcculc-mctal-particlc geometry af- 



fects the EM enhancement and molecule damping rates, 
but also how these parameters in turn influence the spec- 
trum of light emitted by the molecule. In particular, the 
model allows us to simultaneously quantify both scat- 
tering processes (Raman and Rayleigh) and fluorescence 
near metal surfaces, a field that has attracted a growing 
interest in recent years 0, ^| . By applying the model 
to the case of a highly fluorescent molecule situated be- 
tween silver nanoparticles, we obtain an effective Raman 
cross-section of the same order of magnitude as in recent 
smg le-molecule SERS data0,0. In addition, we study 
the effects of a strong incident field, which drives the 
molecule out of thermal equilibrium, and predict that it 
is possible to observe effects such as anti-Stokes Raman 
scattering even at low temperatures. 

Figure ^ schematically shows the main ingredients of 
the model. A molecule is placed on the symmetry axis 
between two spherical, metallic (Ag) nanoparticles. This 
system is illuminated from the side (9 — 90°) by a laser 
with light polarized along the symmetry axis. The scat- 
tered and fluorescent light is collected by a detector also 
placed on the side of the nanoparticle system. 

We treat the molecule as an electronic two-level sys- 
tem with an excitation energy E e — E g = fiVL ge . The 
model also includes one symmetric molecular vibration 
mode with frequency f2 V ib and reduced mass \i. The elec- 
tronic and vibrational degrees of freedom are coupled by 
a Franck-Condon mechanism. As shown in Fig. ^ the 
equilibrium position of the vibrational coordinate is dis- 
placed a distance xq upon electronic excitation. The di- 
mensionless parameter a = xq / y^K]{p^h^) character- 
izes the strength of the electron- vibration coupling which 
ultimately makes Raman processes possible. 

The spectrum (differential cross-section per unit pho- 
ton energy fiuj and solid angle fl) of the light sent out 
by the molecule in the direction of 9 = 90° as a result of 
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FIG. 1: Schematic illustration of (a) the model geometry and 
(b) the oscillator potential governing the vibrational motion 
of the molecule in the electronic ground and excited states. 



both scattering and fluorescence processes can be calcu- 
lated from[nl 
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Here I[ n is the incident laser intensity and M(cu) is the 
EM enhancement factor. The normal-ordered correlation 
function {p^(0)p^(t)) of the molecule dipole moment 
p can be evaluated by means of the quantum regression 
t,heorem|llj once we know the molecule's time-averaged 
density matrix and its equation of motion. In this way 
light emitted both as a result of molecular transitions 
(fluorescence) and as a result of the molecule's oscillating 
dipole moment (scattering) is accounted for. We calcu- 
late the density matrix p describing the molecule dynam- 
ics keeping a finite number iV v ib (usually 4) of vibrational 
levels per electronic level. The equation of motion reads 

i dp/dt = (l/h) [H mol +H', p}+ Ctrp + C ph p, (2) 

where H 1 describes the molecule-laser interaction, £ tI 
and £ p h account for various damping processes which 
we specify below, and H mo i is the molecule Hamiltonian, 
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The molecule-electric-field interaction is a central in- 
gredient in the dynamics that we account for by a —er-E 
term in the Hamiltonian. The electric field E has contri- 
butions both from the incident laser beam and the vac- 
uum fluctuations. The laser field 



El = zEo cos fi^t = zEq [e 
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yields a transition matrix element (in the rotating-wave 
approximation) between two molecule states in the elec- 
tronic ground and excited states, respectively, 

(e; m\H'\g; n) = -M(n L )poE Q e~ mLt f{n, m)/2. (5) 



Here hflL is the photon energy and p — e^jip is 
the transition dipole moment between the electronic 
states. f(n,m) is a Franck-Condon factor, i.e. the over- 
lap f(n,m) — (0;n\xo;m) between state |n) in the 
undisplaced oscillator potential and state |m) in the 
displaced potential, which depends on the parameter 
a = xq/ '2h I '(/iilvib) introduced above [l^. 

For a molecule in free space the interaction with the 
EM vacuum fluctuations yield a decay rate from state 
|e;m) to \g;n) due to spontaneous emission 
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Near the metallic particles the decay rate is modified, 
r 5 „, em — ► |M<i(a;)| 2 r ffn!em , where w = fl ge + (m - n)f2 vib . 
\Md\ 2 = P/Phcc is the ratio of the power emitted by a 
dipole placed at the position of the molecule with and 
without the metal particles present. Usually \Md\ 2 is of 
a similar order of magnitude as |M| 2 , yet the two fac- 
tors may differ substantially because \M<]\ 2 accounts for 
radiation in all directions as well as energy dissipation 
in the metal particles. We calculate M and us- 
ing extended Mie theory The optical properties of 
the particles are represented by a tabulated, local dielec- 
tric function For small (<20-30 A) molecule-particle 
separations d, there are important corrections (~ 1/rf 4 ) 
to the damping-rate enhancement | A/^ | 2 as a result of 
electron-hole pair creation in the particles. To capture 
this we calculate, in the non-retarded limit while apply- 
ing a long-wavelength cutoff, the power P e h dissipated 
by the dipole when placed between two flat Ag sam- 
ples (at the same distances as the spheres) whose optical 
properties are described by a non-local dielectric func- 
tion based on d-parameter theory [16|, and add this to 
the output power Pmi c found in the Mie calculation, i.e. 

\M d \ 2 - (P M ic + Peh)/Pfree. 

Decay and dephasing rates enter the last two terms of 
Eq. Standard quantum optics methods [Tlj yield 
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in the low-temperature limit, (a^j denotes a matrix with 
the only non-zero element kj equal to 1.) is the to- 
tal decay rate from state j to k. It includes the radiative 
and non-radiative processes discussed above as well as vi- 
brational damping due to transitions with a phenomeno- 
logical rate 7 v jb to the nearest, lower level within the 
same electronic state. We also introduce a phenomeno- 
logical dephasing rate 7 P h that enters the last term of Eq. 
(J2J, CphPkj = —il-phPkj provided the electron states of k 
and j differ. Let us stress that 7 P h is brought into the 
model in order to broaden the fluorescence resonances of 
the molecule. In reality an organic molecule has many 
vibration modes, and therefore an almost continuous flu- 
orescence spectrum. Dephasing gives us a broadened flu- 
orescence spectrum even though the model molecule has 
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FIG. 2: Combined Raman and fluorescence spectra for 
a model molecule, with parameter values 7if2 9e =2.3 eV, 
M7 vib =0.1 eV, 7 ph = 10 14 s~\ 7 vib =2xl0 12 s _ \ q=0.5, and 
^di P =l A, placed between two Ag particles with i?=400 A for 
three different molecule-particle separations, and hQ,L=2.45 
eV. Inset: The absorption and Stokes Raman scattering cross- 
sections a a and <t_r as a function of the incident laser photon 
energy for the model molecule in free space. The experimental 
a a for an R6G molecule is shown as a comparison. 
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FIG. 3: Calculated cross-sections as a function of molecule- 
metal-particle distance. The Raman (peaks) and fluorescence 
results have been extracted from the areas marked R and 
F, respectively, in Fig. [21 multiplied by the angular average 
87r/3. The "estimate" curve gives the free molecule Raman 
cross-section multiplied by the relevant enhancement factor 
\M{Q L )\ 2 \M{n L - f2 vib )| 2 . The "total" curve is the total 
cross-section integrated over the energy range 1.5-2.8 eV. 



only one vibrational mode. It has limited impact on the 
resonant Raman scattering as long as 7 P h is smaller than, 
or comparable to, the laser detuning. 

The inset of Fig. [21 shows the absorption cross-section 
a a , and Raman profile (cjr as a function of incident pho- 
ton energy) calculated using the Fermi golden rule for 
the model molecule in free space jl3|. The experimental 
absorption cross-section for a Rhodamine 6G (R6G) dye 
molecule is shown in the same diagram. We have set the 
parameter values cited in the caption to obtain a simi- 
lar spectrum. hfl ge gives the overall peak position and 
?i.n v ib the vibrational quantum, and the value for 7 v ;b is 
reasonable for a molecule at a metal surface, a has been 
chosen to reproduce the shoulder of the spectrum, and 
jph and idip were set to reproduce the width and height of 
the R6G spectrum, respectively. The free-molecule Ra- 
man cross-section is 7 to 8 orders of magnitude smaller 
than a a , and its maximum is blue-shifted compared with 
the absorption spectrum due to quantum-mechanical in- 
terference between processes with different intermediate 
vibrational states for the molecule. 

The main panel in Fig.[21shows spectra calculated with 
the model molecule placed between two silver spheres 
with radius R — 400 A and three different, symmetric 
[i.e. di = rf 2 in Fig. ^ molecule-particle separations. 
All three spectra have a broad fluorescence peak around 
Ml ge — 2.3 eV, which shifts to slightly lower energy for 
d=5 A since the maximum of M is redshifted as the EM 
coupling between the Ag particles increases. In addition 
a number of sharp peaks emerges, due to either Rayleigh 
scattering off the molecule (at 2.45 eV) or Raman scatter- 



ing (red-shifted by multiples of hfl v fo=0.1 eV from 2.45 
eV). ur varies rapidly when the geometry is changed; for 
d = 5 A the Raman peaks are comparable in height to 
the fluorescence background, whereas for d = 15 A the 
Raman peak is barely discernible, cf. Ref. [17J. 

The EM enhancement \M\ grows with decreasing d 
yielding a rapid growth of the Raman signal which in- 
volves both an absorption and emission event and thus 
scales as ~ |M| 4 . Fluorescence also involves both pho- 
ton absorption and emission, but the cross-section in this 
case only scales as ~ |M | 4 /|Md| 2 . The fluorescence in- 
tensity is proportional to the EM enhancement in emis- 
sion, |M| 2 , multiplied by the probability of finding the 
molecule in the excited state. This probability is rela- 
tively insensitive to the enhancement because it is set by 
the ratio (~ \M/Md\ 2 ) between the laser excitation rate 
and the deexcitation rate due to spontaneous emission. 

Figure shows the dependence of the Raman (<7r), 
fluorescence (erp), and total (or) cross-sections on the 
distance d, exhibiting the same general tendencies as 
discussed above, op has been calculated from the area 
marked F in Fig. [21 For the Raman scattering we have 
plotted two curves: one is obtained from the area marked 
R in Fig. [21 but this calculation only works for relatively 
small d, so we also estimate an by multiplying the free 
molecule Raman cross-section by the enhancement fac- 
tor \M(n L )\ 2 \M(n L - n vih )\ 2 . The estimate agrees well 
with the "peaks" result for d ai 15-25 A and, of course, 
gives the true result for larger d. For small d the two 
results differ, and here the peaks result is the true one; 
it includes effects on the Raman scattering of the strong 
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FIG. 4: Stokes and anti-Stokes Raman scattering power (for 
low temperature) as well as their ratio as a function of Ji n 
calculated with the parameter values given in Fig. [5] 

energy dissipation that are not included in the estimate. 

At large d the cross-sections approach those of a free 
molecule. The influence of the particles causes a weak 
interference phenomenon; there is a minimum in (jr and 
op at d «3000 A. For smaller d, both <tr and up in- 
crease as a result of EM enhancement. Over a range of 
distances the enhancement factors roughly scale as 1/d 
meaning that up ~ 1/d 2 and ur ~ 1/d 4 . At d < 30 A, 
more complicated behavior sets in. Resonant enhance- 
ment with resonance frequencies that shift with chang- 
ing geometry occurs, and damping effects become impor- 
tant. These first affect o p which in spite of an increasing 
enhancement M levels off around d = 10 A and even- 
tually decreases, because the molecule is usually deex- 
cited through a non-radiative process. For the smallest 
distances becomes larger than ap, but it eventually 
decreases due to the strong dissipation which damps also 
the coherent oscillations of the molecular dipole moment. 

In Fig. 0] we show results for the power of both Stokes 
and anti-Stokes Raman scattering as a function of in- 
cident laser intensity J; n . For low intensity, the Stokes 
power is linear in I m (constant apt) while the anti- 
Stokes power grows quadratically with I lri . The anti- 
Stokes signal occurs because in an intense laser field the 
molecule can be found in an excited vibrational level 
of the electronic ground state once the rate of elec- 
tronic excitation and deexcitation becomes comparable 
to the vibrational damping rate 7 v it>- The probability 
for the molecule being vibrationally excited is roughly 
C\M(QL)\ 2 aA4>ph/jvib, (C < 1 is a numerical factor), 
and this ratio is « 0.6C with \M(Q L )\ 2 « 1.4 x 10 5 
and ii n =0.5 mW///m 2 corresponding to a photon flux 
^p/i = 10 23 photons/ (cm 2 s). Note that in this model 
the excited vibrational state is pumped mainly by re- 
peated absorption and deexcitation rather than by Ra- 



man scattering 0, For higher intensities both the 
Stokes and anti-Stokes signals saturate and eventually 
decrease. The molecule is driven so hard that its polariz- 
ability becomes time-dependent. This happens when the 
effective Rabi frequency f2_R = M^i^PqEq/Tl becomes 
comparable to other relevant frequency scales, in our case 
the dephasing rate 7 P h- For the parameter values used 
here fl R « 4 x 10 13 s _1 at J in = 0.5mW//im 2 . 

In summary, we have presented a model calculation 
that treats surface-enhanced Raman scattering and flu- 
orescence on an equal footing. We found that, with re- 
alistic parameter values, a resonant Raman cross-section 
of ~ 10 -14 cm 2 can be reached with an EM enhancement 
by 10 orders of magnitude. We also found that for an 
incident laser intensity of ~ lmW//im 2 it is possible to 
get a considerable anti-Stokes Raman signal. 
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